Effective attraction between like-charged colloids in a 2D plasma 
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The existence of attractions between like-charged colloids immersed in ionic solution have been 
discovered in recent experiments. This phenomenon contradicts the predictions of DLVO theory 
and indicates a failure of mean field theory. We study a toy model based on a two dimensional 
one-component plasma, which is exactly soluble at one particular coupling constant. We show that 
colloidal interaction results from a competition between ion-ion repulsion and longer ranged ion-void 
attraction. 
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I. INTRODUCTION 

Recent experiments show convincingly the existence 
of attractions between like-charged colloids immersed in 
ionic solution, in particular in the vicinity of a glass wall 
or when the colloids are confined between glass walls E10 . 
This remarkable counter-intuitive phenomenon is incon- 
sistent with the well established theory of Derjaguin, 
Landau, Verwey and Overbeek (DLVO) and has gener- 
ated various theoretical interpretations. 

It was shown || that counterion correlation forces, be- 
cause of the long range of the electrostatic potentials, are 
usually not pairwise additive. Many-body effects Q and 
Coulomb depletion forces Q have been invoked to explain 
attractive forces between like-charged macroions. An ex- 
act demonstration || proved that the non-linear Poisson- 
Boltzmann (PB) mean-field equation can not give attrac- 
tions in the case of Dirichlet boundary conditions. The 
argument was extended to a broader class of models |7]] 
and boundary conditions ||. Recently, there was a sig- 
nificant new attempt based on non-equilibrium hydrody- 
namic effects 0j. 

In this paper, we solve a 2D one-component plasma 
model exactly at a certain coupling constant following 
the method of Jancovici |ic[| . We have found a semi- 
analytic method which extends the exact solution to the 
case in which colloidal particles are present. We treat 
the colloids in the plasma as empty regions (voids) from 
which the ions are excluded. The colloid may have fixed 
charge at its center. The analytic results obtained for 
small-sized colloids show that while the ion-ion repul- 
sion, which goes like e -fl (R is the separation distance) 
in this particular 2D model, is strong at short distances, 
the ion-void attraction, which goes like R 2 e~ R , becomes 
dominant at large distances. The latter property is not 
normally found in conventional mean field approxima- 
tions, and must come from strong correlation effects. For 
large-sized colloids, we also have numerical results which 
qualitatively agree with the analytic ones. The present 
state of understanding of how to treat reliably fluctua- 
tions beyond mean field theory is not very good despite 



the fact that this is an old problem |ldj-|13j . Thus, our ex- 
act results may serve as a testing ground for more sophis- 
ticated approximation schemes. In addition, since this is 
a 2D model, the results are of interest in connection with 
the problem of interactions among linear polyelectrolyte 
molecules Q. 

In the sections that follow, we shall first review the 
plasma model in section [fl| Then in section III and IV 
we calculate the ion-ion and void- void interactions in the 
plasma system respectively. Section [y] puts the above re- 
sults together by treating the colloid as an empty region 
plus a fixed ion charge at the center. We have relegated 
the technical details as much as possible to the appen- 
dices. 



II. THE PLASMA SYSTEM 

We start by reviewing the exact solution for the one- 
component 2D plasma system Our approach is to 
map the classical plasma system into a quantum Hall sys- 
tem. Such an analogy allows one to describe the interact- 
ing ions in terms of non-interacting electrons occupying 
certain quantum obitals, and subsequently use already 
developed tools. The advantage of this approach will 
become more clear later when we extend the same anal- 
ogy to more complicated systems with fixed ion charges 
or /and empty regions present. 



A. The Quantum Orbital Language 

The one-component 2D plasma system is composed of 
n simple ions embedded in a disk of uniform background 
charge. Each ion carries a charge of —1; the disk has a 
radius r and an areal charge density of po- Throughout 
this paper, we shall express all lengths in the unit of 
In = * — Thus, a unit circle contains unit amount of 

background charge, and the charge neutrality condition 
can be written simply as 



(2.1) 
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In two dimensions, it is convenient to use a complex num- 
ber z — x + iy (in units of £q) to represent the location 
of (x,y). Using Zj to denote the position of the ith ion, 
we can solve the 2D Poisson equation and find the total 
Coulomb potential energy of the system: 
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^ n n 



(2.2) 
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Here the first term describes the attractions between the 
ions and the background charge, and the second term 
represents the mutual repulsions among the ions them- 
selves. 

At inverse temperature (3 — 2, the Boltzmann factor 
reads: 
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Written in the above form, the Boltzmann factor can be 
identified as the square of the determinant of a matrix 
M 



-pv _ 



det Ml 



(2.4) 



with the matrix elements being the lowest Landau level 
(LLL) wavefunctions studied in connection with the 
quantum Hall effect pal: 



Mi 



-M 2 /2 



(2.5) 



The above wavefunction describes an electron occupy- 
ing the ith angular momentum orbital in the LLL. In 
this quantum language, det M is a Slater determinant 
and represents the Laughlin wavefunction at filling fact 
v = 1, describing a fully filled Landau level. Thus, the 
complicated Coulomb interactions are replaced by a sim- 
ple Slater determinant. The determinant representation 
of the Boltzmann factor was first noted by Jancovici |ll]] 
prior to the discovery of the quantum Hall effect. We re- 
capitulate (and then extend) his argument in more mod- 
ern language. 



B. The Partition Function 

Armed with this powerful quantum analogy, we pro- 
ceed to calculate the partition function of the system, by 
integrating out all the ionic degrees of freedom: 



z=n/d 
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(2.6) 



This is nothing more than the norm of the Laughlin wave- 
function, and hence can be computed exactly. Substitut- 
ing eq.(2.4) in cq.(2.6) and expanding the determinant, 
we have: 



n » 

Z = \\ /dz 4 |detM| 2 

i— 1 

n „ n 

=n/ d ^EE(- i ) p+Q n^^)^^)' ( 2 - ? ) 

i=l J {V}{Q} i=l 

where V and Q are permutations of {1 : n} [1 : n is a 
short hand notation for 1, 2, . . . , n] that come fro m th e 
expansions of det M and det M respectively. In eq. ( |2.7| ) , 
we have both summations and integrations which, re- 
spectively represent the Coulomb interactions and the 
thermal averages. It is useful at this point to postpone 
treating the Coulomb interactions until after the ther- 
mal averages are done. Thus, we s witch the order of the 
summations and integrations in eq.(2.7) so that: 



n „ 

^ = EE(- 1 ) P+Q I1 J dz^ Pi (zi)^ Q M) 



{V}{Q} 

n 

{V}{Q} i=l 



(2.8) 



where is the Dirac notation for the inner product of 
the ith and jth LLL orbitals. Notice that the summations 
over {Q} can be identified as a matrix determinant: 

n n 

J2(-i) v+Q tl<n\Qi) = E(-!) K = det °< 

{Q} i=l {7?.} i=l 

(2.9) 

where 1Z = Q o ~p~ l is the composite permutation of Q 
and V , and the overlap matrix O defined by 



Oij = 



M3) 



(2.10) 



is a matrix composed of LLL wavefunction inner prod- 
ucts. Here, just as previously, the complicated Coulomb 
interactions are again replaced by a simp le ma trix de- 
terminant. Finally, we substitute eq.(|2.9|) into ( |2.8| ) to 
obtain: 



J^detO 

{V} 



i\ det O. 



(2.11) 



The immaterial factor of n! will be dropped henceforth. 
Through the quantum analogy, the statistical problem 
is reduced to evaluating the overlap matrix O and its 
determinant. This is particularly simple in the present 
case, because of orthonormality we know that matrix O 
is an identity matrix. 



C. The Plasma Edge 

Near the plasma edge, physical properties are quite 
different from those in the bulk. This motivates us to 
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study the boundary effects separately. For simplicity, 
we are going to be mainly interested in a soft boundary 
condition, which restricts the range of the ions' angular 
momenta rather than their positions. The ions can go 
anywhere in the complex plane, nevertheless, they are 
most likely to be found inside a disk with radius r = \fn 
because their angular momenta are bounded by n [Recall 
that the Zth LLL state has angular momentum of (I — 1) 
and peaks in a shell \Jl — 1 < r < Vl]. The advantage 
of using this type of boundary condition is apparent: we 
can continue to use the same wavefunctions to describe 
the ions near the edge. What is implicitly assumed here 
is that the uniform background charge does not cease 
at the plasma edge; it extends to infinity. We also note 
that in the thermodynamic limit (r — > oo), the circular 
plasma edge becomes locally flat. 

One may also consider putting some 'surface' charge on 
the plasma boundary. But, as indicated below, this effect 
can be absorbed in shifting the position of the boundary. 
We assume the surface charge, if any, does not fluctu- 
ate. Depending on its sign, the boundary charge either 
increases or decreases the number of available angular 
momentum channels. This is more or less equivalent to 
varying the radius r of the plasma disk in a system with 
neutral boundary. Hence, without loss of generality, we 
may assume boundaries are neutral. 

Another type of plasma system we want to address 
in this paper is a plasma strip with two parallel edges 
separated by h. This can be realized in the thermody- 
namic limit by a system of ions with angular momenta 
ranging from (r — -j) 2 to (r + |) 2 . The double edged 
system thus differs from single edged system only in the 
momenta range. To avoid redundancy, we shall use the 
single edged system as example in our calculations. By 
changing the angular momentum limits, the results can 
be easily transferred to the double edged systems. 



III. FIXED ION CHARGES IN A PLASMA 

In this section, we shall apply the same method to 
a slightly more complicated plasma system, with one or 
two ion charge(s) fixed at certain position(s). The result- 
ing partition function is conventionally called the one- or 
two-body density function. The physical significance is 
two-fold. First, we can learn from these density functions 
about how the ions are distributed and how two ions 
are correlated. More importantly, the fixed ion charges 
can be viewed as the limiting case of small-sized charged 
colloids. Thus, the results obtained in this section will 
provide us with many useful insights to the colloidal in- 
teractions. 



A. The One-body Density Function 

First we consider the one-body density function. Wc 
fix the position z\ = w and integrate out the remaining 
ionic degrees of freedom 



Z 



(i) 



it, „ 

n 

i=2 ^ 
n „ 

Yl /dz.ldetM 



= EE(- 1 ) ns ^ I H^ 1 HlI< p 'ift>' 

{V}{Q} i=2 

The summations over P\ and {Q} then reduce to a matrix 
determinant: 

n 

Pi {Q} i=2 

(3.2) 

with the (n + 1) by {n + 1) matrix OW defined by: 



o« = 



O 

$(w) o 



(3.3) 



Here, O (not to be c onfu sed with O^) is the nxn overlap 
matrix defined in ( |2.10| ); 5* = (ipi, ip2i ■ ■ ■ > ipn) is a row 
vector of LLL wavefunctions representing the fixed ion 
charge, and ^ is its conjugate transpose. Fi nally , the 
remaining summations over {-P» (i=2:n)} m e( l-( [Oj ) j us t 
yield a trivial prefactor: 

= -detO« =-(n-l)!detO«. (3.4) 

m (i =2:„)} 

To evaluate the determinant of matrix O^ 1 -*, we apply 
Schur's theorem (discussed in detail in appendix |b|) and 
find: 

=-detO (1) = det O- det [^(w)O- 1 ^ {w)] lxl 

n 

= E^H^H' ( 3 - 5 ) 



k=l 



where we have used the fact that the overlap matrix O 
is an identity matrix. Using eq.(A5), we convert the 
above series into translation coefficients (TCs) defined 
in eq.flA^): 



n i n 



k=l 



k=l 



Now that the partition function is written in the form 
of a TC sum , we may use the TC sum-rules listed in 
appendix A 2. We consider the following two cases. 
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If the fixed ion is in the bulk of the plasma (w <C r = 
y/n), we may take the limit n — > oo first and use TC 
sum-rule ( A12| ) to obtain: 



Z (i) 



1 

7T 



(3.7) 



This tells that deep in the bulk, the ion density is con- 
stant and equal to the background charge density every- 
where. 

Near the plasma edge, we assume the boundary is lo- 
cated at (r + §) and the fixed ion is located at w = re l6 , 
so that the ion is 4 away f rom the bou nda ry. In this case, 
we can use TC sum-rule (A21) in eq.(3.6) to obtain: 



Z (i) a = _ 
2tt 



1 



, h 

1 M 7! 



(3.8) 



Here, the superscript 's' means single edged system (in 
contrast to 'd', which means double edged system) and 
$(x) is the standard error function. As shown in Fig. [l], 
the ion density indeed drops from its bulk value to zero 
upon crossing the soft plasma boundary (at h = 0) . 



B. The Two-body Density Function 

In calculating the one-body density function, we find 
that a fixed ion charge amounts to an additional row 
and column to the overlap matrix. This result can be 
generalized for systems with multiple fixed ion charges. 
In particular, we consider the two-body density function 
with z\ — w and 22 = w fixed. The result is: 



=-detO( 2 \ 
where O^ 2 ) is an (n + 2) by (n + 2) matrix: 



Q(2) 



O &(w) &(w) 
tf(ttf) 




(3.9) 



(3.10) 



Similarly using Schur's theorem and eq.(A5), we can 
write the correlation function in terms of TC sums: 



Z {2) = det O • det 



*(w) 

-i 2 



O 1 [ *t(to) &(w) ] 



Lfe=i 



E 3 ^ 

k=l 



(3.11) 



In the bulk of the plasma (w -C r), we assume w = 
so tha t the two fi xed io ns are R apart . Plugging TC sum- 
rules ( A12 ) and ( A13 ) into eq.( 3.11 ), and using eq.( A7) 
for Tn(ii?), we obtain: 



^) = i[l-|Tn(ii?)| 2 ] = i(l-e-« 



(3.12) 



The above result tells that the Coulomb repulsion be- 
tween the two ion charges is heavily screened by the other 
ions and becomes short ranged. Its decay has a Gaussian 
form, which is much faster than the exponential form 
that DLVO theory predicts. 

Near the plasma edge, we assume the boundary is lo- 
cated at (r+ -|) and the fixed ions are located at w — re 10 
and w, where 9 = M so that the fixed ions are -| away 
from the e dge a nd R apart from e ach o ther. Using TC 



asymptotically large R): 



sum- rules (A.21) and (A24) in eq.(3.11), we obtain (for 



z (2)s = _±_ 

Air 2 



1 + $ 



V2 



2e 



ttR 2 



(3.13) 



The above result shows that the ion-ion correlation func- 
tion has a power law decay near the plasma edge. This 
is fundamentally different from the bulk behavior, and is 
due to the dipole moments induced by the sharp cutoff 
in angular momentum. 

Finally, in the double edged system where the plasma 
edges are lo cated at (r ± -|), we may use TC sum-rules 
(A2E) and (A31) in eq.(3.11) to obtain (for asymptoti- 
cally large Rj: 



z (2)d = 



1 



$ 2 



h 

V2 



2ttR 2 



sin 2 (hR) 



(3.14) 



Here, we find the correlation is not only long ranged, but 
oscillating with a period of 2irh~ 1 as well. 



C. Divalent Ions 

A divalent ion can be viewed as composed of two simple 
ions occupying the same place. Here we want to calculate 
the two-body density function for a pair of divalent ions 
embedded in the monovalent plasma. Let us assume z\ = 
Z2 = w and 2:3 = Z4 = w are fixed, where w = i% so that 
the two divalent ions are R apart. Naively one may follow 
eq. ( 3.10] ) and write down an (n + 4) by (n + 4) matrix: 



6(4) 



O *t( Zl ) $t( Z2 ) ^t( Z3 ) ^t( Z4 ) 
#(zi) 

*(z 2 ) 

v{z 3 ) 0000 

*(z 4 ) 



(3.15) 

However, the matrix O^ 4 ' defined above is singular be- 
cause z\ = Z2 (or Z3 = 24). This singu lari ty actually 
originates from the Boltzman factor in eq. ( [2.3] ) , where we 
incorrectly include the self interaction between z\ and 22 
(or 23 and 24). To get around this problem, we must ex- 
tract the singular factors from the matrix determinant. 
More precisely, we want to consider the following limits: 
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detO» S lim lim 

z 2 ^zi z 4 ^z 3 z 2 — 21 r Z4 — Z3 r 



(3.16) 



To construct matrix O^ 4 ) from d^ 4 \ we subtract ^(zi) 
in t he (n + l)th row from ^(#2) in the (n + 2)th row in 
eg. ( 3. 15 ). The resulting (n + 2)th row then yields a factor 
of (z 2 - Zi): 

lim [*(* 2 ) - = (z 2 - zi)*'(zi), (3.17) 

22— >Zl 

where VP' = (V4> V^j ■ ■ ■ , i>'n), and V*' is the first derivative 
of ip. More precisely, following the usual procedure for 
the Hilbert space of analytic functions [ fulfill , we take the 
derivative to act only on the analytic part of tp and not on 
the Gaussian factor. The latter represents the interaction 
with the background charge and does not possess the self- 
interaction singularity. 

The factors of (22 — ~zi ) and | z 4 — z 3 1 2 can be extracted 
via similar procedures and we are left with a matrix that 
looks like the following: 



o< 4 > = 



o 

*(«;) 
*'(w) 



&(w) &\w) ^(W) *t'(E7) 



(3.18) 



Applying Schur's theorem and using various TC sum- 
rules listed in appendix A 2, one finds: 



-R 2 /2 







=-detO< 4 > 

1 

1 Re~ R2 / 2 e- fl2 /2 

e- R2 / 2 Re~ R2 / 2 1 R 
e- R2 / 2 R 1 + R 2 
2 



i det 



1 



i? 4 e 



-R* 



(3.19) 



To compare the correlation functions of divale nt ion s and 
simple ions (in the bulk case), we plot eq.( |3.12| ) and ( 3. IS ) 
in Fig. |^. Clearly, the divalent ions has a larger 'exchange 
hole' near its origin. 



IV. VOID COLLOIDS IN A PLASMA 

In this section, we introduce two empty colloidal voids 
into the plasma system. We take the colloidal particle to 
be a circular disk with finite radius a, inside which ions 
are excluded. The presence of the voids clearly destroys 
the system's azimuthal symmetry. As a result, different 
angular momentum orbitals intersect, and the overlap 
matrix is generally not diagonal. It is the off-diagonal 
matrix elements that make the colloidal interactions non- 
trivial. 



We shall fix th e two voids at w and w, respectively. Fol- 
lowing eq.( 2.11 ), we can write down the partition func- 
tion of the system as the determinant of the overlap ma- 
trix: 



Z = detO. 



(4.1) 



Because the regions occupied by colloids are not accessi- 
ble to the ions, their contribution to the ove rlap m atrix 
should be excluded. Thus, in contrast to eq.( 2.10 ), here 
we have: 



Oi 



~ (i\j)w - (i\jh 



(4.2) 



where (i\j) w means a partial inner product integrated 
ove r a circle wit h radius a centered at w. According to 
eq.(A6), eq.(4.2) can be explicitly evaluated as [ p~7| : 



Oi 



j2"/i T *i( w ) T Ji( i 



1=1 



(4.3) 



where 7; is the incomplete gamma function defined in 
eq.Q. 



A. Analytic Results 

The off- diago nal matrix elements make the determi- 
nant in eq.( |4.l| ) difficult to analyse. However, if the col- 
loidal siz e a is small, we can treat them perturbatively. 
From eq. ([A3]) we know 7; ~ 0(a 2t ) when a — > 0. Thus, 
to the lowest order in a, we can truncate the infinite sum- 



mation in eq.(4.3) at I = 1: 

0^ « 8 l3 - 7l \T %1 {w)T 3l {w) + Ti^wfTj^w)] . (4.4) 



If we define a 2 x n matrix J as 

Tu(w) Tzi(w) ■■■ T nl (w) 
Tu(w) T 2 x(w) ■■■ T nl (w) 



eq^^) may be rewritten in a matrix form: 
O = 1- jt 7l J. 



(4.5) 



(4.6) 



According to the corollary of Schur's theorem (see ap- 
pendix pi eq.(|B3|)), we have: 



det(l-jt 7l j) nxn = det(l-J 7l jt; 



2x2 ' 



(4.7) 



The determinant of a 2 x 2 matrix is easy to calculate 
and the result is (dropping terms with order higher than 
7i): 

Z = detO«l-2^ \T k i{w)\ 2 ^j 7 i = 1 - 2 7l , (4.8) 



where we have used the TC sum- rule (A12). Notice that 
the above (first order) result is i?-independcnt. This is 
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no accident, because in order for the colloids to know 
the separation distance R, the Green's function needs to 
be integrated over both colloidal regions, which yeilds a 
term proportional to a 4 , appearing in at least the second 
order in perturbation theory. 

In the second order approximation, we must keep terms 
proportional to 71, 7 2 and 72. The procedure is quite 
similar to the above, and the result is: 



|T fcl H| 2 71-2 VlIkHI 2 72 




T kl {w)Y 



k=l 



(4.9) 



Compared to eq.(|3.6|) and ( |3.11 ), we find that the first 
and second order terms are respectively proportional to 
the one- and two-body density functions. Hence the col- 
loidal void-void interaction does not differ very much 
from the ion-ion interaction we studied earlier in Section 

in. 

In the p lasm a bulk , we assume w = i-j 
sum-rules (A12) and (A13) in eq.(4.9) to obtain 



I-271-272+I 



,-R* 



7 2 ; 



(4.10) 



near the plasma edge(s), we assume w = re (9 — , 
and express eq.(|4.9|) in terms of TC sum-rules results: 



Z>>* a 1 - 2^71 " 2 ^22 d 72 + {A\?f - \B s {i 



(4.11) 



where the explicit forms for A s ' d and B s ' d for a symp toti- 
cally large R can be found in eq.( A21 ), ( A24 ), ( |A28| ) and 



B. Numerical Results 

For large-sized colloids, the analytic results based on 
the small a-expansion may not be appropriate. We re- 
sort to numerical methods and use LU factorization with 
partial pivoting to evaluate the matrix determinants. We 
find that (essentially) exact results can be obtained from 
relatively small matrix sizes. 

For colloidal voids immersed in the bulk of the plasma, 
our numerical results are presented in Fig. |^. For a vari- 
ety of sizes we studied, no attraction is found, confirming 
our analytic predictions. 

Near the plasma edge(s), we find it more convenient to 
use square shaped colloids in numerical calculation. We 
assume the squares are of size 2a, centered at re and 
re -10 respectively. [The plasma edge(s) are located at 
(r ± 4).] De finin g 5 = we can calculate the overlap 
matrix in eq.(4.2) as following: 



Oij — &ij 



r+a r-O+S p i+]-l e -p 2 e iU-i)4> 

dp I dej) , h cc. 



(4.12) 



In the limit i ~ j ' ^ r 2 — > 00, the above integral can be 
evaluated asymptotically. The result is 



Oi 



where, r± 



c(o,6) 



r(^) 

Vmm 



[$(r+)-$(r_)], (4.13) 



_ (r±a) 2 -i±2- 



-■ C(0 5) = cos [('~J') 6 '] sin[(i-j)S] 
' v ' / (i— jW 



-2 j — (%-}) 

(C(9, 5) should be interpreted as — in case of i — j). In 
practice, we keep increasing r until the numerical result 
becomes independent of r, which is found to occur ap- 
proximately at r = 60. The results shown in Fig. |] (for 
single edged system) and Fig. || (for double edged system) 
are computed using r = 100. In both figures, when the 
colloidal voids are far away from the edge(s) (h = 6.0), 
the curves are similar to those in Fig. |^, namely, the re- 
pulsive interaction has a Gaussian decay. As the colloids 
become closer to the plasma edge(s) (smaller h), the re- 
pulsion starts to develop a long tail, in agreement with 
our analytic results. 



V. CHARGED COLLOIDS IN A PLASMA 

We now consider charged colloids in the plasma. In the 
simplest case, the charged colloid is taken to be an empty 
disk plus a fixed ion charge at the center. The ion-ion in- 
teractio n a nd void-void interaction have been studied in 
Section III and IV respectively; both are repulsive. How- 
ever, as we shall see in this section, the ion-void interac- 
tion is usually attractive, and has a longer range. The 
competition between the repulsion and attraction results 
in a richer behav ior for the charg ed colloidal interactions. 

Following eq. (|3.9| ) and ( 3.10| ), we have the partition 
function as follows: 



= 



-detO (2) 
O *t(u,) 
*(iu) 









(5.1) 
(5.2) 



where, VP = [tpi, ip2, ■ ■ ■ , ip n ) is a row vector representing 
the fixed ion charge at the colloi dal c enter; the n x n 
overlap matrix O, according to eq.(4.3), is non-diagonal: 



E H T ti ( w ) T ji ( w ) + c - c - 
.1=1 



(5.3) 
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A. Analytic Results 

If the colloidal size a is small, we can obtain an ana- 
lytic expression. One may truncate the infinite series in 
eq.(0) at I = 1: 



O, 



<% - 7i \T iX {w)T jX {w) + T iX (w)T jX (w)] . (5.4) 



However we notice th at t he (n + l)th and (n + 2)th row 
of matrix O^ 2 ) in eq.(5.2) can be written as 



O 



(2) 



o 



(2) 

(n+2)i 



J* 
T jx (w) 



(5.5) 
(5.6) 



If we multiply (|J) by ^K^ x T ix (w), multiply fl5.6| ) by 
V / '' : 7i^ 1 ji(' u; ); an d add them t oget her to the iih row, the 
off-diagonal part of Oij in eq. ( |5.4| ) is cancelled competely. 
Thus the final answer does not contain 71. This is true 
even if we do not make any truncations in eq.(5.4). Phys- 
ically, it originates from the fact that the first angular 
momentum channel is occupied by the fixed ion, so it is 
not accessible to the free ions. Another way of saying 
this is that the fixed ion serves as a 'Laughlin quasi hole' 
in the quantum electron system. 

To obtain a non-trivial result, we consider the second 
angular momentum channel 1 = 2: 



O 



m ~ 5ij - 72 \T i2 (w)T j2 (w) + T i2 (w)T j2 (w)] . (5.7) 



Applying Schur's theorem, we have: 



Z (2) = det O • det 
= Z\ ■ Z 2 . 



*(w) 



cr 1 [ &(w) &(w) 



(5.8) 



We first concentrate on the second factor Z 2 . It describes 
the two fixed ions (represented by wavefunctions W and 
^) interacting through the plasma medium (represented 
by matrix O -1 ). For small 72, the matrix O 1 can be 
approximated by: 



o: 



Sij + 72 [T l2 {w)T ]2 {w) + T l2 {w)T j2 {w)} , (5.9) 



and Z 2 can be calculated to be (suppressing the (w) ar- 
guments) : 



E Ir« 



vfc=l 



+272 




Ei T «i 



E TkiTk 2 



k=l 



E T 'i E T ^T k 2 E T ^ 

\k=l ) \fc=l / \k=l ) 

/ n \ / n \ / n ^ 

E t m E T « r « E T * T « 



vfe=l 



\fc=l 



\fe=l 



(5.10) 



We can interprete the above result by thinking of the col- 
loidal voids as a perturbation to the plasma syste m plus 
fixed ions. The first square bracket in eq.(5.10) is the 
result for the unperturbed system, namely, the ion-ion 
interaction via the plasma medium. It is p recise ly the 
two-body density function we obtained in eq.(3.11). The 
second square bracket in eq.(5.10) is the first order cor- 
rection, namely, a fixed ion on one site interacts with the 
colloidal void on the other site. This ion-void interaction 
is usually found to be attractive, and since it involves 
higher order TC sums, its range is longer than that of 
the ion-ion repulsion (see results below). So far, we have 
missed another first order effect, namely, the fixed ion 
interacts with the colloidal void on the same site. The 
excess background charge underneath the void should re- 
duce the effective ionic charge somehow. It turns out 
that this renormalization effect is captured precisely by 
the first factor Z\ in eq.(5.8): 



Z x =dctO« l-2^|T fc2 | 2 72. 



(5.11) 



fe=i 



In the bulk of plasma, we assume w = '\% and use TC 



sum-rules (A.12) and (A13) to obtain 



'(2) 



(1 - 2 72 ) 1 - c- H + 2R 2 e~ H 72. (5.12) 



In the above equation, the first term is the renormal- 
ized ion-ion repulsion and the second term is the ion- void 
attraction. The attraction is proportional to R 2 e~ R , 
which has a longer range than the repulsion. This fea- 
ture is not found in mean field approximations, thus must 
be caused by the strong correlation effects. Notice that 
the void-void repulsion and the renormalization of the 
ion-void attraction are both of second order, so they do 
not appear in eq.(5.12). The repulsion and attraction re- 
spectively dominates the short distance and long distance 
behavior, and the minimum free energy 0F = —\nZ^ 
occurs where the first derivative vanishes: 



dR 



= 



R* = J— ■ (5.13) 



Fig. H plots the calculated free energy vs. R for vari- 
ous wducs of a. Note that when 72 is small (a = 0.6), 
eq. fl5.13| ) predicts that the location of the free energy min- 
imum is far away from origin, in which case the attraction 
is strongly suppressed by the decaying factor e~ R , thus 
becoming virtually invisible. 



Near t he p lasma edge, we use TC sum-rules ( A21 ) 
through (A.25) to obtain (for asymptotically large R): 



7 



z (2)s 



2irR 2 



+272 



+ ^1 



(Ai 1 h* + A5 3 )e- h ' "^e~^ 



'2-k 



2irR 2 



/1 ll /1 22 



(5.14) 



Simil arly for do uble edged system, we use TC sum-rules 
(|A2SD through (|A||) to obtain 



7(2)d 




■ sin 2 (hR) 



2irR 2 

[A^h 2 cos 2 (hR) + A d 22 sm 2 (hR)}e- h ' 



2vri? 2 



«) 2 ^2 2 }, 



(5.15) 



B. Numerical Results 

Our numerical results for charged colloids in the bulk 
case are plotted in Fig. ^[ When a is small (a — 0.6 — 0.8) 
the numerical results agree very well with the analytic 
ones presented in Fig. |6| However, when a is large (a = 
1.0), the two differ significantly in that the numerical 
free energy has an additional secondary maximum. This 
peak reflects the even longer ranged, second order effects 
which are not accounted for in the analytic studies. 

We also tried putting double ion charges at the center 
of each colloids. This case is very simila r to the divalent 
ions which are studied in Section Iff C . Our numerical 
results are presented in Fig. ||, where the colloidal sizes 
are chosen to be \[2 times as large as those in Fig. |?], 
so that the two figures may be comparable. We see that 
the qualitative behavior is very similar. The attraction 
minimum stays further out. 

Near the plasma edge, we use squares to represent col- 
loidal particles as we did in Section IV B . The numerical 
results for two different sizes of colloids are presented in 
Fig. I and Fig. |h]. In case I (a = 0.8, Fig. §), when the 
colloids are far away from the edge (h = 5.0) they repel; 
as they become closer to the edge (h = 2.3) they start 
to attract each other within a certain range (R < 2.25), 
beyond that range they still repel; at a particular point 
(h = 1.5), the range of the attraction seems to extend 
to infinity; finally, when they are too close to the edge 
(h = 0.8), due to insufficient screening, a long range re- 
pulsion is found. In case II (a = 0.9, Fig. [l(]), the colloids 
are larger in size and attract each other even in the bulk 
(h = 5.0); the attraction becomes deeper (h — 2.0) and 
longer ranged (h = 1.3) when the colloids move towards 
the edge. 



Fig. |TT| shows the free energy of colloids in a double 
edged system. When the two boundaries are widely sep- 
arated (h — 5.0), the colloids repel; as the separation 
decreases, the free energy develops a long oscillating tail. 
These are consistent with the analytic result s. No te that 
the oscillation 'wavelength', according to eq.( [5.15 ), is set 
by the strip width h, and is independent of the size a of 
the colloids. This conclusion is confirmed by Fig. [l2| 

Finally, Fig. [l3| shows the results for doubly charged 
colloids near a single plasma edge. They are qualitatively 
similar to those shown in Fig. Il0|. 



C. Where Does Colloid Stay? 

To better understand the charged colloidal interac- 
tions, especially near the plasma edge, it is instructive 
to learn where they want to stay in the plasma. We put 
a charged colloid at w = re ie , which is | away from the 
plasma edge at (r + 4), and calculate the partition func- 
tion to be: 



z (l)s = 

o« = 



det (1) 

O &(w) 
V(w) 



where the n x n overlap matrix O is: 



Oij — Sij 



^27iT u (w)Tji(w). 



(5.16) 
(5.17) 



(5.18) 



It can be shown that the first angular momentum channel 
is blocked by the fixed ion charge so that the final answer 
does not depend on 71. For small a, we truncate the 
infinite summation in eq.(5.18) at I = 2 and use Schur's 
theorem to find: 





n 




n 


n 


2 




l-El T ^| 2 72 




i=l 


'y] Ti\Ti2 
i=i 


72 



= A° n 



2tt 



-A S U A S 22 72 . 



(5.19) 



Fig. [IT] shows the numerical results done for a square 
charged colloid. We plot the free energy (3F as a function 
of h for three different sized squares. In all cases, we 
find the colloid is attracted to a free energy minimum 
near the plasma edge (at h = 0). For h < 0, the free 
energy quickly blows up. This tells that the colloid does 
not want to flee away from the plasma edge, because its 
fixed charge is attracted by the background charge and 
wants to be in the lowest available angular momentum 
channel. The optimal distance h* is in the range in which 
attraction between colloids occurs. 
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VI. SUMMARY 

In this paper, we considered a strongly coupled {(3 = 2) 
2D one-component plasma model and calculated the ef- 
fective interaction between two colloids both in the bulk 
and near boundaries. We found that due to ion-void in- 
teraction, the colloids may attract despite the fact that 
ion-ion and void-void interactions are both repulsive. 
Near the plasma edge, the attraction effects become more 
prominent. Though our model is oversimplified, the ex- 
act results provide useful insights into the problem, and 
can serve as a testing ground for more sophisticated ap- 
proximate treatments. 
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be expressed as a linear combination of these basis wave- 
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APPENDIX A: PROPERTIES OF LOWEST 
LANDAU LEVEL WAVEFUNCTIONS 

Some important properties of lowest Landau level 
(LLL) wavefunctions are summarized here. They are fre- 
quently used throughout the paper. 



1. Orthonormality and Completeness 



The LLL wavefunctions defined in eq.( |2.5| ) are or- 
thonormal: 



= / dzipi(z)ipj(z) = 5 i: 



(Al) 



Sometimes, we want to calculate a partial inner product 
(the overlap integrated over part of the entire plane), and 
the result is generally not diagonal. However, if the inte- 
gration region has azimuthal symmetry, the conservation 
of angular momentum guarantees an exception. As an 
example, let us take the region to be a circle with radius 
a centered at origin (denoted by the subscript), 



WJ)o 



P d P / d ^ /T,-\Tf\ = 6 * jlU 

Jo 7r v /r( i )r(j) 



where the imcomplete gamma function 7, is defined as 
the following: 



2. Translation Coefficients 



The linear coefficients T^s defined in eq.(A4) are 
ter med translation coefficients (TCs). Setting z = in 
eq.(A4), we find: 



T»i (to) 



(A5) 



Thus the LLL basis functions are just special cases of 
TCs. Another useful application of TC is to calculate the 
partial inner product on a circle with radius a centered 
at w ^ (denoted by the subscript): 

{i\j)w= / dzip i (z + w)i()j{z + w) 

J\z\<a 

N <Q k,i 



= Y J T lk (w){k\l) a T 3l { 

k,l 
00 

= ^2jiT u (w)T j i(w). 



(A6) 



1=1 



The TCs can be evaluated explicitly. We multiply 
ipk(z) to both sides of eq.(A4), and integrate z over the 
entire plane. On the right hand side, this picks out one 
particular from the summation because of the or- 
thonormality: 

T lk {w) = jdze^-^^iiz + w) 
e -M 2 /2 



e -M 2 /2 



Ky/T(i)T(k) 



dze -\z\ -^z k ~ 1 (z + 



dze-l^ 2 



w 



9 



£ 



k* i— k 



i(i,k) 



£ 



(-i)>i 



e 2 



i+fc-2s 



(-+f) 2 _ 

■Kj lim £ T fel (re ie )^(re ie ) 

fc=l 
(>'+f) 2 

lim V T fe4 (re i9 )T fc ,(re ie ). 



5f 



^ r(s)(i-s)!(fe-s)! 



(A7) 



(A16) 



(A17) 



fe=i 



Where, in the 3rd equation we used Taylor and binomial 
expansi ons, and in th e 4th equation we used orthonor- 
mality jAl). Eq.(]A7l) immediately tells that the TCs 
have the following inversion symmetry: 



Tij{w) = Tji{-w) = Tji(-w). 



(A8) 



Another important property of TC arises when one 
composes two successive translations: 



Tij (wi + w 2 ) = c 



J2 T ik(wi)T kj (w 2 ). (A9) 



fc=i 



Here, the superscript 's' denotes single edged system and 
9 = tP-. The difference between the bulk TC sum and 

2r 

edge TC sum is as following: treated as a perturba- 
tion, colloids in the bulk mainly perturb the low an- 
gular momentum channels. Such a perturbation has a 
discrete characteristic, namely only a small number of 
channels are strongly affected. Colloids near the bound- 
ary, on the other hand, perturb the high angular mo- 
mentum channels. The perturbation has a continuum 
nature, namely there are a large number of channels 
(r 2 — hr < k < r 2 + hr) being weakly affected. As a 
result, the evaluation of edge TC sum inevitably requires 
one to convert the discrete summation into a continuum 
integral. For example, let us calculate Af x : 



Or using eq.(A8), this becomes: 



T ij (w 1 +w 2 ) = ei^' w ^J2T kt (~ Wl )T k3 (w 2 ) (A10) 



4 s 



if 



k=l 

oo 



lim 



lim 



£ 

k=\ 



„2fc-2 



W-LW-2—W 



lW2) Y, T ki{-^) T kM (All) 



r(fc) 

r'+hr dx 



2x In r-\-x — x In x—r 



(A18) 



k=l 



As a sim ple application, we set — Wi 
eq.(|ATo|): 



where we have defined x = k — 1 and used Stirling's 
asymptotic formula: 



W 2 



T{x + l) » V2^e- X+Xlnx . 



(A19) 



lim y T kt {w)T k3 (w) = Tn (0) = , (A12) 



k=l 



or set wi = W2 



The integal in eq.( A18 ) can be evaluated using the saddle 
point approximation, in which one expands the exponent 
near its minimum at x — r 2 up to the second order: 



in 



in eq.(All) 



A\± « lim 



r +hr 



dx 



lim y T) 



a | — | Tkj 



k=l 



Tij(iR). (A13) substituting y 



o \/2nr 
— , one finds 



-e 2 



Note that if we interprete Ty as the matrix elements of 
T, then eq.(A12) has a simple matrix form: 



Ji lx — 



dy _ 



y A /2 _ 



1 



2tt 



1 + $ 



V2 



(A20) 



(A21) 



T\w)T{w) = 1, 



(A14) 



showing that ma trix T is unitary. In the same matrix 
notation, eq. ([A8|) is rewritten as 



Where $(a;) is the standard error function. Other TC 
sum rules can be obtained through quite similar proce- 
dures. We omit the tedious derivations and simply list 
the results below: 



Tt(io) = T(-w) 



(A15) 



A s — - 



which along with eq.( A14 ) makes sense since T(— w) 
should be t he inve rse o f T(w). 

Eq.(|AT2]) and ( |A13| ) are examples of TC sum-rules. 
They are frequently used in analysing the bulk of the 
plasma. Near the edge, we often encounter TC summa- 
tions of a different type: 



.4 



12 



B{ x = 



Biry — 



1 + $ 

e- h2 / 2 

e -h'/2 



h 

71 



he~ h2 ' 2 



i(hR-%) 



2nR 



he~ h2 / 2 . 



i(hR-f) 



2vri? 



(A22) 

(A23) 
(A24) 
(A25) 
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Note that the above results are correct only asymptoti- 
cally for R oo (R=2rd). 

In a double edged system, we encounter the same type 
of TC summation, except that the lower limit is changed 
from k = 1 to k = (r — |) 2 : 

(r+|) 2 

4=lim V T kl (re w )T kj (re ie ) (A26) 

J r — >oo * — » 

B% ee lim g ^(re ie )T fc ,(re ie ). (A27) 

fc=(r--|) 2 

Correspondingly, we have: (the superscript 'd' denotes 
double edges): 



^2 = * 
id 



ft 

71 



fte"'* 2 / 2 



A? 2 = 0, 

#11 



r>d 
B 12 



/2ttR 
he~ h2 / 2 



2nR 



V2n ' 

sm{hR)e i( - hR -i\ 
cos(hR)e i{hR -^. 



(A28) 

(A29) 
(A30) 
(A31) 

(A32) 



APPENDIX B: SCHUR'S THEOREM 

For any four matrices A rxr , B rxs , C sxr and D sxs , if 
A is not singular, Schur's theorem states 



det 



A B 
C D 



detA-det(D-CA" 1 B) 



(Bl) 



The proof is straightforward. We consider the following 
matrix identity: 



l r x r 






A 


B 




CA 1 


Isxs 




c 


D 





A B 

sxr D — CA _1 B 



(B2) 



Taking the determinant of both sides of eq.(B2), eq.(Bl) 
follows immediately. 

Schur's theorem has a very useful corollary: for any 
matrices 3 r xs and K sxr , one has 

dct(l - JK) rxr = det(l - KJ) SXS . (B3) 



We can prove this via very similar method. First we 
consider 



1, 



J sxr i sxs 



K 1. 



(1 JK) rxr rxs 
K l, x . 



(B4) 



which tells 



det 



L r xr 

K 1 



j 

sxs 



det(l — JK) rxr ; 



on the other hand, if we consider 



0,. 



i rxr u rxs 

-K l sxs 



we obtain 
det 



L r X r • 

K 1. 



± r xr u 

sxr l-KJ 



L r X r • 

K 1. 



det(l - KJ) 



(B5) 



(B6) 



(B7) 



A combination of cq.(|B5|) and (B7) results in cq.(p33j). 
Letting J = 71K''' in eq.(p33|), we obtain eq.(4.7). 
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by 

/3F = -IndetO 

= -trln{l - [T(m)GT t (w) +T(ra)GT t (W)]} 

OO 

= ^ Ur[T{w)GT t (w) + T^GT^tt)]* 



fc=i 



]T-trG fc + FW 

fc=i 

2/3F + 0F(w) 



where, (3Fo = — lndet(l — G) represents the free energy 
of each individual void, and F(w) can be viewed as the 
interaction potential. 




FIG. 1. The one-body density function near the plasma 
edge. 
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FIG. 2. The comparison between the two-body density 
functions for divalent ions (the solid line) and simple ions 
(the dash line) in the bulk. 
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FIG. 3. Numerical results for the free energy of colloidal 
voids of radius a in the bulk of plasma separated by a distance 
R. 




FIG. 4. Numerical results for the free energy of empty col- 
loids of radius a = 1.0 and located a distance h/2 away from 
a single plasma edge as a function of their separation R. 




FIG. 5. Numerical results for the free energy of empty col- 
loids of radius a = 1.0 separated by a distance R in a plasma 
strip of width h. 
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FIG. 6. Analytic results for the free energy vs. separation 
distance R of singly charged colloids in the bulk of the plasma. 



FIG. 9. Numerical results for the free energy vs. separation 
distance R of singly charged colloids near a single plasma 
edge; case I: a = 0.8. 
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FIG. 7. Numerical results for the free energy vs. separation 
distance R of singly charged colloids of radius a in the bulk 
of the plasma. 
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FIG. 10. Numerical results for the free energy vs. separa- 
tion distance R of singly charged colloids near a single plasma 
edge; case II: a = 0.9. 
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FIG. 8. Numerical results for the free energy vs. separation 
distance R of doubly charged colloids of radius a in the bulk 
of the plasma. 
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FIG. 11. Numerical results for the free energy vs. separa- 
tion distance R of singly charged colloids near double plasma 
edges; a — 0.8. 
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FIG. 12. Numerical results for the free energy vs. separa- 
tion distance R of singly charged colloids near double plasma 
edges; h — 1.5. 
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FIG. 13. Numerical results for the free energy vs. separa- 
tion distance R of doubly charged colloids near a single plasma 
edge; a = 1.2. 
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FIG. 14. Numerical results for the free energy for one singly 
charged colloid near a single plasma edge. 
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